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ANNALS OF M ATHEMATICS. 

Vol. III. June, 1887. No. 3. 

COPLANAR MOTION OF TWO PLANETS, ONE HAVING A ZERO MASS. 
By Dr. G. W. Hill, Washington, D. C. 

The supposition that two planets circulate about their central body in the 
same plane enables us to dispense with two differential equations of the second 
order in the general problem of three bodies. The further supposition, that the 
mass of one of them is too insignificant to have any sensible effect on the motion 
of the other, enables us to consider the motion of the latter as known and as 
taking place according to the laws of Kepler. Hence, in this case, the two co-or- 
dinates of the planet of zero mass are the only unknowns ; and they are given by 
two differential equations of the second order. These suppositions have, approxi- 
mately, place in several cases in the solar system, but I have more especially in 
view the motion of the satellite Hyperion as disturbed by the action of Titan. 
My object in this paper is simply to point out a method of proceeding, which 
may, I think, be advantageously employed in this case. 

Employing the usual notation x, y, r, for the rectangular co-ordinates and 
radius vector of the planet whose motion is to be determined, xf , y' , r', for the 
corresponding quantities belonging to the acting planet, m' the mass of the latter, 
and M the mass of the central body, the differential equations of motion will be 

d 2 x_3Q_ dy_9Q 
dt 2 ~ 3x' dt 2 ~ 3y • 

where ii, the potential function, has the following expression : — 

n- M 4. „/ T I xhc^/f\ 

~ v \x 2 +f) + M Li/[(* - x'f + {y -y'f] r» J ■ 

The co-ordinates of m' satisfy the differential equations 

d 2 x' , M+m' , d 2 y' , M+m' . 

We can, without any loss of generality, assume that the axis of x is directed to- 
ward the lower apsis of m'. Then the integrals of the last-stated differential 

equations are 

x' = a' (cos e' — e'), y' = a' j/(i — e' 2 ) sin e', 



66 HILL. COPLANAR MOTION OF TWO PLANETS, 

where e' is derived from the equation 

n't + c' = e' — e' sin s', 

a', ^', c' being constants, and »' being the equivalent of A 7s — 

It is desirable to know what the differential equations determining x and y 
become when expressed in terms of any other variables. For this end Lagrange's 
canonical form of the equations serves very conveniently. Let the new variables 
be u and s, and employ the subscript Q to denote the complete differential co- 
efficient with respect to / of any variable to which it is attached. Then T stand- 
ing for ^ (xi + y*) expressed in terms of u, s, u lt s lt Lagrange's canonical form 
of the equations is 

d9T _2T _3ii d 3T_2T __2ii 
dt 2u t 3u 3u ' dt 3sj 3s 2s 

. , 3x , 2x 2x 

As we have - r i = *r a i + T" J i + a: » 

3u 3s 2t 

,, — 9 ? „ _u dy c j- $ y 



we get 

¥ L U«J + U«J V 1 + l9u 3s + 3u 2. 



, i r r^v-f { 3 a is* + [ 9x$x + ^ ^ 

+ * L UJ + l&J J 1 + U«^ ?* & 

, (3x 3x ,3y2y) f f 2x^\ 2 , fc>^ H 

It is very plain from the form of Lagrange's equations that if the variables 

u and .f were so assumed that one of them, u for instance, should disappear at 

once from the expressions for T and ii, we should have an integral of the prob- 

d 2T 2T 

leva. For then -.- «■- = O; and, integrating,, = a constant. This selection, in 
dt2u l 2th 

a theoretical sense, is always posssible, and in as many essentially distinct ways 
as there are first integrals of the problem, which, in the present case, are four. 
But, although it is easy, in innumerable ways, to make U depend on one variable, 
it is not so eas;/ to make the six factors of the general expression for T depend 
solely on the same variable. And, when we inquire what equations must be 
satisfied for this, we find that they are essentially the same as those which are 
satisfied by the Eulerian multipliers. Hence, nothing is gained by approaching 
the problem from this side. 

I propose to take u and s so that 

x — fix' it -f- f>y's, y = f>y'u — fix's. 
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where p denotes a function of t supposed known, but, for the present, left inde- 
terminate. From these equations may be derived 

r 2 = p 2 r n (u 2 + -s 2 ), x'x -f- _/j» = />r' 2 w. 

Hence the potential function, in terms of u and s, becomes 

o_ i_ f M m^__ , 2 1 

U ~ pr> L ,/ (« 2 + ^) + j7 {(« - P - ') 2 4- *] m '° U \ ' 

In the general expression for T we substitute the values 



Sx , Sy . dx . Sy , 

Tu = F x > Tu = W' Ts=&> Ts = -> ox > 

3x_d( I >x') d{p/) 3y_d{p/) d(px') 

St ~ dt + dt • St ~ dt dt 

The result is 

T= y 2 r' 2 (it 2 + V) - * n »' i/(i - O P 2 («*i - «i) 

For the sake of brevity we may write, /Zj, /z 2 , /z 3 being known functions of /, 

T= \h x (u* + s 2 ) - k 2 (us, - su,) + \h % (u 2 + s 2 ). 

This, substituted in Lagrange's canonical form of the differential equations, gives 
as the equations of the problem, 

dt V h dt J + 2 ' h dt *** + dt S ~ Su ' 

d f , ds'} , du dh, , Sii 

dt K dt) 2 dt dt 3 Ss 

Let us now adopt a more general independent variable than the time. Call- 
ing this £, let dt= dd^, in which d may be regarded as a function either of/ or £. 
The second supposition will be the more advantageous. In either case as we 
obtain, on integrating, u and .$ as functions of f , it will be necessary to have the 
values of f which correspond to given values of the time, and thus the inverse 
function will have to be considered. Then, in terms of the new independent va- 
riable, 

d [>hdi±\ , 2& d L _ d M, 9J0S) 

d'z U < J + 2 < 3 < Su ' 

d (I h ds^ du dk 2 dhu _2{0Q) 

dt U < J "~~ 2 dt ~ dC "** ~ ST ■ 

We can now consider how p and 6 should be assumed in order that the dif- 
ferential equations may be most simplified. In the first place it appears impor- 
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tant that the potential function Q should be freed from the independent variable 
£. This is accomplished by putting p = i. In the second place it seems we 
cannot readily do better than take the eccentric anomaly e' of the attracting 
planet as the independent variable £. Then 





r' r' 
dt = —7—. de' and 6 = — r - 7 
a'n' a'n 


\lso we have 




h x \d = a n n'{i 


— e' cos e'). K = a n n V( i — e'~), 6 




~ = ^- f , M+m'=a' i n' 2 . 



e'% Oh 3 = a' 2 n'(i +e'coss'), 



pr 



For the sake of simplicity let the signification of Q be changed, and, putting 

m! , 

-jr-f-. 7 = v, let 



,/(** + /) ^ >/[(*_ i)* + ^] 
Then our differential equations take the form 



^[(.-,'cose')g] + 2l /(i-^)§-(i+,'coseO« = ^, 

*. [ (I _ e > cos «0 §] - 2i/d ~ *' 2 ) J, - (I + • cos «') , = ^ . 

It will be noticed that the potential function Q is, by this assumption of va- 
riables, completely freed from co-ordinates expressing the position of the attract- 
ing planet ; and that the two factors i — e' cos s' and I -f- e' cos e', very simple 
functions of the independent variable e', are the only evidences of the position of 
this body in the differential equations. And, of the four elements of its orbit, e' 
is the only one we have to deal with. 

We propose now to see whether the introduction of elliptic co-ordinates will 
bring about any simplification in the problem. Supposing 



. 2 



let — hr H rr= l > and — rr H rT- =I ' 

be the equations of a confocal ellipse and hyperbola, a x and # 2 being constants, 
and l x and 1, the new variables destined to take the place of u and s. By elimi- 
nating x£ from these equations we obtain 

a, — <?, , 
— x \ = I ; 



K + *i) K + 4s) 
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whence j| >, + A) K + ^ 1 

\ L a x — a 2 J 

The expression of x 2 in terms of X x and A, is obtained from this by simply inter- 
changing a x and a 2 . Thus 

x = |r (gi±A)(gL±Ai "| 

2 \ L «2 — a, J ' 

We now proceed to find what i? becomes in terms of l x and ^ 2 . By taking 
the sum of the squares of the last two equations we get 

x \ + x i = a \ + a 2 + A + 4 
Thus far a x and a 2 have been left indeterminate, but we now assume 

a 2 — a x = \- 
Then u 2 + s* = (x 2 + |) 2 + *, 2 

= 2a 2 -f- Aj + l 2 + 2 T /[(« 2 + ^) (« 2 + A 2 )] 

= [i/(« 2 + ^i) + v / K + ^)] 2 , 

^/O 2 + /) = ,/(«, + X x ) + j/(« 2 + 4), 

(„_,)' + ** = (*,-*)* + *,» 

= 2« 2 + a, + ,* 2 — 2|/[(« 2 + ^,) (« 2 + ;, 2 )], 
v [(« - i) 2 + /] = /k + ;,) - !/(« 2 + ;. 2 ), 
« = 2 1 /[(rt 2 + ;,) (« 2 + ; 2 )] + 1. 

For the sake of brevity we will now put 

V{a 2 + *i)=p. V( a 2 + K) = q- 

Then it is plain Q may be written 

n I — v , v 

We have now to deal with T. By taking the logarithms of the values of 
x x and x 2 , and then differentiating, we obtain 

dx x _ dl x ( dl 2 



*\ ' a x + l x a x + ;, 2 

dx^ _ dl x <fl a 

2 .. — T" i)i 



i, a 2 + /[ « 2 + / 2 



2 

I J1 _ 

2 
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Whence may be derived 

On substituting in the factor of dk x dl 2 the values of x? and .r 2 2 it vanishes, and 
the expression takes the form 

4 (dx? + dxi) = , /! ~ k - ■ ■ dl? + . VW^T ;, dh\ 

Or, in terms of/ and <?, we have 



du 2 + ^ 2 = £ — t <*/* + ^ — 4- 2 ^ 2 - 

/ — i ^ — i 

In like manner we get 

uds-sdu = (f, + q)[ y j l^i) Jf-yj (f£^J <*?]• 
The former expression for 7" was 

7-=i(i-^cosO^^-i/(i-0 ' ftfc ^; J<fe + i(i+^cosO(« , + ^); 

hence, if we abbreviate by putting 

7--Ki-^co.o[(i+^)^+[i+i]g] 

- l/(l - O (/ + *) [« ^ - I %\ + i (I + e' cos «') (/ + qf. 
T and Q are somewhat simplified if we adopt variables p and a, such that 

/> + q — p, p — q = °- 

Also, for the sake of brevity, put 



Then we have 



- 1/(1 - n p{k% + hYe\ + i (1 + ^' cos «') p\ 



Q = i V + - — W + W. 
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By this transformation ii is considerably simplified ; but, as more than off- 
setting this, T is rendered complex. As the expression for a in terms of these 
variables is 



tt - jrL-^-^n 



it will be perceived that h and k are trigonometrical functions of the angles of the 
triangle whose sides are I , [>, and <x, which might have been anticipated from geo- 
metrical considerations. Thus it appears no advantage would result from the 
employment of elliptic co-ordinates. 

Returning, therefore, to the quasi-rectangular co-ordinates u and s, it seems 
some advantage would be gained if we adopt a new system of co-ordinates u and 
s, such that the new system is expressed, in terms of the old, as follows : — 

u = u -\- s y ( — i), s = u — S ]/( — i). 

We can also adopt the trigonometrical exponential corresponding to the arc s' 
as the independent variable. Calling this £ = e e ' l '^~ n , an operator D is adopted, 

equivalent to £ -=- , so that D . C = iC- 
- dC - 

In terms of the new variables, ii has the expression 

l/{us) !/[(«— I) (j— i)] 2 v ' 

And the differential equations are 
D {[i -VC+ C- l )\Du} + 2,/(i -e>*)Du + [i + \c> {' + C" 1 )] « = - 2 ~~ , 

d {[i ~yc + r 1 )]^} - 2 1 (i -e' 2 )Ds + [i + 1/(: + r 1 )] *= - 2 C ^ . 

Only one of these equations need be actually employed, as either can be obtained 
from the other by changing the sign of j/( — i). We have 

c ii _ i — v v 

c'/t? _ I — v v 

For the purpose of integrating these equations, we may adopt the method 
of indeterminate coefficients ; and we may employ, as proper to represent the 
values of u and s, the infinite series 

" — — • <* ,, ;, (■ ^ , 

c V a r — '« —./'=' — '• 
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Here i,j, and k denote positive or negative integers, zero included ; and the sum- 
mation must be extended so as to include all values for i, j, or k from — oo 
to + oo . The a and c, c' are constants and functions of the four quantities e', v, 
a and e ; a and e being two of the four arbitrary constants introduced by integra- 
tion. The two remaining arbitrary constants serve only to complete the two 
elementary arguments which belong to the attracted planet, and, in this method 
of integration, they can pass unnoticed. 

If we suppose that the orbit of the attracting planet is circular, the differen- 
tial equations reduce to the very simple form 

(D-lfs = -2^. 

And, in this case, an integral can be found. For multiplying the first by Ds, and 
the second by Du, the sum of the equations, thus multiplied, is an exact deriva- 
tive. Integrating, we get 

Du Ds + us + 2J2 = 2(7, 
C being the arbitrary constant. 

This integral equation may be combined with the differential equations in 
such a way that one of the terms, regarded as the most difficult of expression in 
a developed form, may be eliminated. For example, if this is taken to be the 

term — -=-, *—. ™ of ii, the equations serving to determine the a may be 

!/[(«— i) (s — i)] • h is r 

taken to be 

(s- i)D(D + 2)u + IDhDs+ ±^u + $(u-v)(s-v)+ C=o, 

(« -i)D(D-2)s + \DuDs + y^ 3 s + *(«-*) (*-»)+ C= o, 

in which the constant C is not identical with the former C. One of these equa- 
tions suffices, as the other is a consequence of it. The difference of these equa- 
tions is simpler than either of them, and may be of use. It is 

D [(« - i) Ds - (s - i) Du - 2 (m - i) (s - i)] = ^t^ 3 (« — *)• 

In attempting to derive periodic series for the co-ordinates of Hyperion, it 
appears to me that it will be easier, in the first instance, to assume that Titan 
describes a circular orbit. And, in the next place, to assume that the perturba- 
tions are periodic functions of the mean elongation of the two bodies. And, as 
it may very easily happen that the terms, depending on the second and higher 
powers of the disturbing force, may quite alter the values of the coefficients, it 
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will be well to employ the method of mechanical quadratures. Starting Hype- 
rion from its line of conjunction with Titan, and at right angles to this line, with 
an assumed velocity, trace out its path until the elongation, between the two 
bodies, amounts to 180 . Then, if Hyperion is again moving at right angles to 
its radius vector, the velocity at the start has been rightly assumed. But if not, 
one makes another trial ; and, by interpolating between the two results, a velocity 
is obtained which will more nearly bring about this condition. And continued 
repetition of these trials will enable us to discover, with all desired approximation, 
the velocity which fulfils this condition. When the path of Hyperion, corres- 
ponding to this velocity, has been traced out, it will be easy, by the well-known 
processes of mechanical quadratures, to assign the periodic series representing 
the co-ordinates of the satellite under the supposed conditions. 

When this is done, corrections to the co-ordinates, proportional to the first 
power of the satellite's proper eccentricity, can be obtained by the integration 
of a linear differential equation. By comparison of these with observation an 
approximate value of this proper eccentricity will be obtained ; a thing to be de- 
sired as we seem to know next to nothing about it at present. Also one will be 
enabled to decide whether the motion of the mean anomaly is more rapid than 
that of the mean longitude, as has been asserted, without sufficient reason as it 
seems to me. 

As illustrating this point, suppose that our moon, instead of having an eccen- 
tricity about 0.055, na< ^ one about o.ooi. Then the variation would be the pre- 
vailing inequality, and the moon would appear to be in perigee always about 
syzygies, and in apogee about quadratures. In consequence the perigee would 
appear to retrograde with reference to the sun as fast as the moon advances with 
reference to the same body. And yet the relation between the motion of the 
argument, denominated the mean anomaly, and the motion of the mean longi- 
tude, would be nearly the same as it is at present. But the position of the peri- 
saturnium of Hyperion has been concluded from its observed shortest and longest 
radii vectores. This is allowable only when the inequality, called the equation 
of the centre, is the overpowering one. 

After the terms, proportional to the first power of the eccentricity, have been 
obtained, those factored by the second, third, etc., powers, can be derived by in- 
tegrating differential equations of the same character. 

In applying the process of mechanical quadratures to the motion of Hype- 
rion, one will meet the difficulty of the uncertain value of the mass of Titan. 
But this cannot be avoided ; an assumption must be made, and the results after- 
wards corrected by comparison with observation. 



